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ABSTRACT 

We discuss the Hamiltonian formulation of the Schwinger proper-time method of calculating 
Green functions in gauge theories. Instead of calculating Feynman diagrams, we solve the cor- 
responding Dyson-Schwinger equations. We express the solutions in terms of vertex operators, 
which consist of the zero modes of the vertex operators appearing in string theories. We show 
how color decomposition arises in this formalism at tree level. At the one-loop level, we arrive 
at expressions similar to those obtained in the background gauge. In both cases, no special 
gauge-fixing procedure is needed. 
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It has always been difficult to compare theory with experimental data of strong interactions, 
because calculations of higher-order diagrams in non-abelian gauge theories are very compli- 
cated. In recent years, considerable progress has been made in simplifying these calculations, 
by introducing new techniques (spinor helicity method color decomposition [Q, supersym- 
metry identities [§], recursion relations 0], electric circuit analogy [§, etc). An interesting 
approach has been developed by Bern, et ai, and is based on string theory They exploited 
the organization of string theory diagrams to formulate rules for the calculation of diagrams 
in non-abelian gauge theories. Although these rules are inspired by string theory, they can be 
formulated in pure field theoretical terms. The question then arises whether similar organizing 
principles can be derived from field theory. This is an important question, because such princi- 
ples will take us beyond the one-loop order in perturbation theory. They may also permit the 
inclusion of fermion lines in diagrams, which is not an easy task in the present string theoretical 
formulation. 

A more direct approach to the re-organization of Feynman diagrams in non-abelian gauge 
theories has been taken by Strassler [0]. He showed that one-particle irreducible functions in an 
external gauge field may be written as quantum mechanical amplitudes of particles moving on 
world lines. To this end, he employed the Schwinger proper-time formalism |^ and expressed 
the results in terms of one- dimensional path integrals. This approach seems to rely heavily on 
the existence of a Lagrangian which is either quadratic in certain fields, or can be written in 
this form. Thus, it appears to be limited to one-loop diagrams. f\ 

To go beyond effective actions, we shall use the Schwinger proper-time formalism in order 
to solve Dyson-Schwinger equations perturbatively. The solutions cannot always be written as 
path integrals of particle Lagrangians with a simple interaction term. For explicit calculations, 
it seems to be advantageous to use the Hamiltonian formalism, and express the amplitudes as 
expectation values of a string of vertex operators, which are reminiscent of the old operator 
formalism in string theory. 

We start by considering the case of an external gauge field. In this case, the solution to the 

t Extensions to this approach have also been studied H. 
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Dyson- Schwinger equation exponentiates, which explains why it is possible to express it in terms 
of a one-dimensional path integral. We shall discuss this exponentiation in the Hamiltonian 
formalism in order to generalize it later. For definiteness, consider a scalar field (j) coupled to 
an external gauge field = A'^°'T°', where are the generators of the algebra of the gauge 
group. The classical equation of motion for is 

{td^ + gA^{x)y(f){x) = 0. (1) 

In quantum field theory, the propagator A^^{x,x') = satisfies the corresponding 

Dyson- Schwinger equation, 

{td^ + gA^{x) fA{x, x') = 6\x - x') , (2) 

where we have suppressed group theory indices. In order to simplify the notation, we shall 
consider the abelian case first. We may solve Eq. (j^) with the aid of the Schwinger proper time 
T as follows. Define operators and satisfying the standard commutation relations, 

[p/^ , x''] = ir^^"" . (3) 

Let them describe a system whose evolution is governed by the Hamiltonian 

H={p, + gA,{x)f . (4) 

If |g) and \q') are eigenstates of the momentum operator, p^, the function 

t/(g,g';r) = (g'|e-^ng), (5) 

which is a matrix element of the evolution operator, U{t) = e"*^"^, satisfies the equation, 

tdrU{q, q'- r) = {q, + gA,{x)fU{q, q'; r) , (6) 

where = —id/dq^. Integrating over the proper time r, and using the orthogonality of 
momentum eigenfunctions, {q\q') = {2ny6'^{q + q') (treating both momenta as outgoing), we 
obtain 

[q, + gA^{x)f / drf/(g, g'; r) = -i{2Tif5\q + q') . (7) 

•J 



Comparison with Eq. (0) shows that i dTU{q, q'; r) is the scalar propagator in momentum 
space. So, instead of using Feynman diagrams to calculate the propagator, we may instead 
calculate the evolution function, f/(g, g';r). For certain gauge fields, this function can be 
calculated exactly. However, we want to calculate U (g, q'; r) for general fields, because we are 
interested in dynamical gauge bosons. We shall therefore use perturbation theory. 

By splitting the Hamiltonian as H = Hq + gHj, where 

Ho = p^p^ , Hj = {p^ , A^^} + gA^A^ , (8) 
we can write the evolution operator in terms of time-ordered products as 

U{t) = f^tl9}l r dt,--- r dt^T{Hj{t,)---Hj{t^))=Te-''fo'^'"'^) , (9) 

where time evolution is now governed by Hq. Thus, x^{t) and p^(t) satisfy the Ehrenfest 
relations 

^ = [i/o, x1 = 2tp^ , ^ = [i/o, p1 = . (10) 
The solution to these equations is 

x^'{t) = 2ip>'t + , p''{t) = pf" , (11) 

where x'^ and p'^ are the position and momentum operators, respectively, at time t = 0. We 
can therefore write 

f/(g,g';r) = (g';0|re^/o^'^*^^(*)|g;r), (12) 

where the state \q;t) satisfies the Schrodinger equation, idt\q;t) = Ho\q]t). It is clear from 
Eq. ([121) that U {q, q'; r) can be expressed as a quantum mechanical path integral, with ap- 
propriate boundary conditions at times t = and t = t. This is possible because the terms 
involving the interaction Hamiltonian, Hj, have neatly exponentiated. We wish to generalize 
this to the case where such an exponentiation is not possible, so we shall use the operator 
formalism to calculate the various matrix elements in Eq. (|12]) . By expressing (t) in terms of 
x^(r) and x^(0), 

x^{t) = ^x'^iT) + il-^)x^iO), (13) 
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we can write Hj in terms of operators T4(0) and Vfe(T), where the vertex operator Vk{t) is 
defined by 

Vk{t) = e**^-"(*) . (14) 

This operator consists of the zero modes of the tachyonic vertex operator in string theory. Ma- 
trix elements containing insertions of this operator can be calculated by using the commutation 
relations 

VkAr)VkM = e-'^'-'^^VkMVkAr) , [^fci(O) , T4.(0)] = [H,(t) , 14,(t)] = , (15) 
and the properties 

{q-0\Vk{0) = {q + k-0\ , Vk{T)\q- r) = \q + k; r) . (16) 

We have now assembled all the tools which are necessary for the calculation of U{q, q'; r). By 
expanding in the coupling constant, 

U(q,q';r)^U^'\q,q';T)+gU^'\q,q';r) + ... , (17) 

we obtain 

U^'\q,q';r) = {q';0\q;r) = ^^'^^27:^5% + q') , (18) 

whose integral over r is the propagator for a scalar, 1/q^. To first order in the coupling constant, 
we have 

U^'\q,q';r) = r dt{q';0\{p^ , A'^mq^r) . (19) 

The projection of this functional onto a plane wave, — e^e^^'^ — e'*T4, can be written in the 
form 

U^\q.(^Mr) = {q-q'),£dt{q';0\Vk{t)\q;T) . (20) 
After some straightforward algebra, we obtain 

Ul'\q,q',k;T) = {2n)'S\q + q' + k){q - q'), J^^ dte-^^'-^'^^'>'' . (21) 

Upon integrating over t and r, we recover the standard expression, -^{q — q')fj,-^- This can 
easily be obtained by amputating the two scalar legs (which amounts to multiplying by q^q'^), 
and integrating by parts twice. 



4 



To second order in the coupling constant, we have 

U^'\q,q';r) = (g';0| dh dh{p, , ^^(^2)}^ , A''(h)} + £ dtA,(t)A^(t)^ \q;T) . 

(22) 

Projecting onto plane waves, we obtain the evolution function 

U^Jiq, q', h, ^2; r) = {q'; 0| (/J" dt, I^' dt^ip,, Vk,{h)}{pu, VkAh)} + {f^ ^ 

(23) 

+ V^.uIodtVk,{t)Vk,{t))\q;r) . 
After some straightforward manipulations, we can bring this into the form 

Uj^Jiq, q', h, t) = {2q + A;i){^,(2?' + ^2).} /J" dh J^' ^^^e-^,^ 

(24) 

which can be shown to agree with the expression obtained from Feynman diagrams. To obtain 
the one-loop scalar contribution to the self-energy of gauge bosons, we need to sew the two 
scalar legs together, after amputating one of them. We therefore need to calculate the quantity 

U^XkiM.r) 0^,q\2nrS\q + q')trUj^J{q,q',h,h;r) . (25) 

To calculate Uj^{ki,k2;T), first we integrate by parts (using g^e""^*^ = ~^^~^^^)' obtain 
U^Jih,h;T) = -/^(2g + A;0{,,(2g + M.}/o^^^e-^'-^-*)^^-^'^+'^^ 

(26) 

^'llJ.v J (27r)4 ^ 1 

where we discarded a total derivative with respect to r, and omitted a (5-function imposing 
conservation of momentum ((27r)^5^(/i;i -|- k2)). It is now easy to integrate over the momentum 
by completing the square in the exponent. The result is 

U^J {k^, k,- r) = { dt (^77,. + ^(1 - 2^^rk^,k,.) e-'=?*(i-*/-) - ir?,.| . (27) 

By integrating the first term by parts, we can write Uj^ {ki, k^-, r) in a manifestly gauge invariant 
form, 

U^Xki. k2; t) = -^(^M^i. - ri,X) [ dt{l - 2\fe-^l'^^-'l^^ . (28) 
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Next, we consider the case of a fermion, ijj, for which the propagator in momentum space, 
A(g, q') satisfies the Dyson-Schwinger equation 

U+g4{x))A{q,q') = (27:)' S\q + q') , (29) 

where x'^ = —id/dq^. We can express A{q,q') in terms of A'{q,q'), where 

A{q,q')=^A'{q,q'). (30) 

Then A'(g, q') satisfies the equation 

{q' + g 4(x) A'{q, q') = {2nY5\q + q') . (31) 

To calculate A'{q,q'), we may follow the same procedure as before, except that now the inter- 
action Hamiltonian is simpler, 

Hj=4{x)^. (32) 

Notice that we did not have to use the second-order formalism for fermions, which would amount 
to introducing the propagator A"{q,q'), defined by 

A{q,q') = U+g4{x))A"{q,q'), (33) 

instead of A'(g, q') (Eq. (|30|)). The latter leads to more elegant expressions, but the interaction 
Hamiltonian is more complicated. As a result, our computations will be simpler than those 
in the path integral approach for fermions [0. Of course, the final expressions in the two 
approaches agree with each other. 

The nth-order contribution to the evolution function U{q, q'; r) can be written as a time- 
ordered matrix element, as before, 

U^^\q,q'-T)= r dt^--- f"" dt^{q'-M4{tn)^---m)^W.r) . (34) 
JO Jo 

To give an example of an explicit calculation, consider the one-loop fermion contribution to the 
self-energy of gauge bosons, 

f/;;i°°P(A:i, k2;r) = J-0^J (^^'(2vr)^5^(g + q')tTU^^J{q, q\ h, h; r) , (35) 
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where (c/. Eq.(0)) 

Uj^J{q,q',h,k2;T) = £ dt,j\t2{q' -,01 {'Jf,Vk,{t2) ^uVk.ih) ^+ (/i ^ z/)) |g;r) . (36) 
Manipulating this expression as in the scalar case, we obtain 

Ul'Jih, h; t) = ^ih,h, - r^,,kl) dt{l - i^Y-e~^l'^^-'/^\2TiY5\h + h) , (37) 

which is in a manifestly gauge-invariant form without any further manipulations. 

Next, we turn to the calculation of amplitudes involving dynamical gauge bosons. In order 
to calculate amplitudes, we need to fix the gauge. We shall choose the Lorentz-Feynman gauge, 
in which 

d^Al = , (38) 

and the propagator for the gauge field is 

{A;ik)Alik')) = ^5%.l(2vr)V(fc - k') . (39) 

First, consider the gauge coupling of a scalar field. Let 

Tf^{q,q\k) = {<l^%)<l>^iq')A;ik)) (40) 

be the three-point function, with the two scalar legs attached. To lowest order in the coupling 
constant, it satisfies the Dyson-Schwinger equation (in matrix notation) 

qXiQ, q'; k) = gT^iq, , e'^'^jAiq, q') , (41) 

where = —id/dq^ and A{q,q') is the scalar propagator, corresponding to the evolution 
function U^^^q, q';T) (Eq. ([ISD). Eq. ( PD may be solved by introducing the Schwinger proper 
time. To use previous results, notice that U^\q,q' ,k]T) (Eq. (p!9D) obeys the equation 

z9,f/««(g, q', k- r) = q'Uj^'^^iq, q', k; r) + T^{q^ , e^^-^}U^'\q, q'; r) , (42) 

which is a consequence of Eq. (|^) . It follows that 

POO 

r^(g, q'- k) = tg drUj^'^iq, g', k; r) . (43) 
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Therefore, the three-point correlation function corresponds to the first-order evolution function, 
U^^\q, q'] r), even if the gauge boson is a dynamical field. This simple correspondence extends 
to higher orders in the abelian case (as may be seen recursively, by using the appropriate Dyson- 
Schwinger equation for (0(g)0(g')A'^^ ■ ■ ■ A'^")), but fails for non-abelian groups, because of the 
gauge boson self-coupling. Yet, higher-order non-abelian diagrams can still be described by 
proper-time integrals. We shall demonstrate this in the case of a pure Yang-Mills theory. It 
will then be clear how one can extend the formalism to include matter fields. 

Classically, the gauge field satisfies the non-abelian Maxwell equation, which is conveniently 
written in matrix form in the Lorentz-Feynman gauge, 

UA^ = 2g[A, , d^A^] + g[d'^A, , A^] - g[A% d^A,] + g^[A'^ , [A, , A^]] (44) 

To write down the Dyson-Schwinger equation for the three-point coupling, 

Gf:^{K h) = {A';{h)Al{k2)A;{h)) , (45) 

we need to attach the leg with momentum ki. The function 

G;%{ku k2, k,) = "l^T^Gl^iku k,, k,) (46) 

satisfies the equation 

klG;l%{k^,k2,h:) = gT'T%{r^^,{k^-k2), + r^.p{2k2 + k^^ , 

(47) 

where = —id/dki, and 6"''^ri^^A is the propagator for the gauge field (Eq. (P^D). The dots 
denote a term which is obtained from the first one by a permutation of the two amputated 
legs (color decomposition [Q). It can easily be seen that the solution to this equation can be 
written in the form 

G;%{ku k„ k,) = kiA;%{ku k„ k,) + kiA;%{k„ k„ k,) , (48) 

where is the color-ordered partial amplitude, 

A;%{k,,k,,ks) = igT'nlj^ dr dt{k,;0\n%{k2;t)\ks;T) , (49) 
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and 

t) = Vpu{2p,Vk{t) - Vk{t)pp) + Vup{2Vk{t)p^ - PpVk{t)) + Vp^Mt) , Pu} ■ (50) 

Manipulating this expression as before, we recover the standard expression, 

Gf;^(A;i, h) = tg r'^Mh - h)p + Vupih - h)p + Vp^^ih - h),) + ... . (51) 

Similarly, the four-point amplitude with one leg attached satisfies the equation 
klG*Xy^{ki, ...,h)= igj ciW^i-^ { [T''' , T^'] ((2A^'(x)9'^A;^'(x) + 9'^A^'(x)A;^'(x) 

-A''-ix)d,Aiix))At{h)A';,ih)Aiik4)) 

+g[T''' , [T^' , T'^']](A^'-(x)A^'(a:X(x)A^(A;2)A^(fc3)A^(A;4)) } . 

(52) 

The various terms are conveniently organized by color ordering. The result is 
kfGX^ih, = tgT\r]^,{h - k^). + 7],,{2k2 + k{)^ - 7]^^{2ki + k2)u)e''''''Al%i{ki, ks, k^) 

+ige'^-^--A*X{h, k2, k^)T\T^^p{h - k^)^ + 77,^(2/^3 + h), - r]xA2k, + k^)p) 

+g^T^T^T\7],p7]^^ + 7]^p7]xu - 27]^p7]^,)e'^^-^e'^-^-^A{ki, k^) 

+ ... , 

(53) 

where = —id/dki. It is interesting to note that the last term on the right-hand side 
corresponds to the Feynman rule for the four-point gauge coupling obtained by string theoretical 
techniques This is obtained in our formalism without any special gauge fixing procedure. 
The solution to Eq. ( ^3)) can be written in terms of partial amplitudes (c/. Eq. (^9D) which are 
integrals over proper time r, of 

u;^Jpxiku...,k,) = {krM[r, dU S',^ dm%{k2-MW^hk,-M) 



+ IodtnZx{k2,ks))\k,;T) + 



(54) 

^pupx y 

multiplied by the appropriate generators of the gauge group, where T^^i^}p{k]t) is given by 
Eq. (|0D, and 

27]xpVpu)Vk,{t)Vk,{t) ■ (55) 
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Clearly, Eq. ( p^ does not correspond to a quantum- mechanical interaction Hamiltonian. Yet, 
the four-point amplitude can be written in terms of the vertex operator (Eq. (0)). 

Loop diagrams can also be calculated in the same fashion. The one-loop contribution to the 
gauge boson self-energy in a pure Yang-Mills theory, r°^(/ci, ^2), may be calculated by using 
the Dyson- Schwinger equation 

kfrilih, k2) = igj d^x e'"'-^ { [T^' , T'']{(2Ald'' A^l^ + d'^A^A^^ - A^'''d^Ai){x)Al{k2)) 

g[T'' , [T^' , T''']]{A'''^{x)Ai{x)A'l{x)Al{k,)) 
+ [T'\r']lp'\x)id,u^\x)Al{k2))] , 

(56) 

where the last term accounts for the effects of the ghosts. It can be calculated using the 
appropriate Dyson- Schwinger equation. Eq. (^) may be written as 

/c2r;t(fci, k2) =igj dqliqX" -{q + kirv^ne'^' ^Glitiq, ^2) , T^'] 

+ g[T''' , [T''' , T'']]e'''^-''e'''^-''rif,^A{q, q') (57) 

+ e*'=-g^G'fet^)^(g,g',A;2)} , 



where dq = 7|;;^7f^(2vr)^(5'^(g + q'), and = —id/dq^. The solution i 



(27r)4 (27r) 



IS 



roc 

^tuiki, k,) = g^r^'-'f^'-' dTUl;'^°^{k„ k,- t) , (58) 

where 

f/^;'°°P(fci, k,- t)=J ^,q\q: 0\nl-'°°^{k,, k,- T)\q- r) , (59) 

and 

r) = dt, /o*^ dt^(n^^l{k^- h)nf"^(j\k,; t,) + nf\k2; hWf\k,- 1,] 

+ dtr]^M,{t)Vk,{t) , 

(60) 

where Ti'-^^ is given by Eq. (^Of), and 

nlf''\k;t)=p,V,{t) (61) 
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is the interaction Hamiltonian corresponding to the ghosts. After performing the operator 
algebra and integrating over the momentum as in the scalar case, we obtain (c/. Eq. (p?])) 

rf^ih, h) = ^r^'^'f'^'^' /o°° rfr{ dt{j,r^,, + 1(1- 2^fkr,h, 

(62) 

which can be cast into a manifestly gauge-invariant form by integrating by parts (c/. Eq. (|28|)). 
This expression is similar to the one obtained in the path-integral formalism, if the background 
gauge is used 0. However, we arrived at Eq. (^) without any special gauge- fixing procedure. 



Matter fields can also be included in this formalism. The Maxwell equation (|^) is modified 
by the addition of a current, which for a scalar field is 

j^ = 0*r'^(9^-*^A^)0-c. c. . (63) 

This term introduces a correction to the Dyson- Schwinger equation for the gauge boson self 



energy (Eq. (^)). It can be calculated similarly to the other terms. It is also possible to 
derive expressions for graphs with external matter legs, by using the Dyson-Schwinger equation 
corresponding to the Maxwell equation, if the graph contains external gauge bosons, or the 
classical equation of motion for the matter field (Eq. (|lD). We have demonstrated this by 
calculating the three-point coupling (Eq. (^OD). The calculation of higher-order amplitudes can 
be done by a recursive use of Dyson-Schwinger equations. 

In conclusion, we presented a method for the calculation of Green functions in non-abelian 
gauge theories, which is based on the Schwinger proper time formalism , rather than Feynman 
diagrams. Using Dyson-Schwinger equations, we arrived at expressions which resemble those 
derived by string theoretical techniques 0. Similar expressions for one-particle irreducible 
functions have also been derived by expressing the effective action as a path integral for a single 
particle moving on a world-line . The method presented here is more generally applicable than 
the one based on the path integral. It can be applied to tree diagrams as well as loops. Higher- 
order diagrams can be expressed in terms of lower-order ones. Thus recursively, any diagram 
can be written in terms of vertex operators which consist of the zero modes of the vertex 
operators in string theory. At tree level, color decomposition arises without any special gauge- 
fixing procedure. At the one-loop level, we obtained expressions similar to those obtained in the 
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background gauge. It is hoped that this formahsm will lead to more compact and manageable 
expressions for higher-order corrections to Green functions. 
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